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ABSTRACT 

We present a general framework to treat the evolution of one-point proba- 
bility distribution function (PDF) for cosmic density 5 and velocity-divergence 
fields 9. In particular, we derive an evolution equation for the one-point PDFs 
and consider the stochastic nature associated with these quantities. Under the 
local approximation that the evolution of cosmic fluid fields can be character- 
ized by the Lagrangian local dynamics with finite degrees of freedom, evolution 
equation for PDFs becomes a closed form and consistent formal solutions are 
constructed. Adopting this local approximation, we explicitly evaluate the one- 
point PDFs P(S) and P(&) from the spherical and the ellipsoidal collapse models 
as the representative Lagrangian local dynamics. In a Gaussian initial condition, 
while the local density PDF from the ellipsoidal model almost coincides with the 
that of the spherical model, differences between spherical and ellipsoidal collapse 
model are found in the velocity-divergence PDF. These behaviors have also been 
confirmed from the perturbative analysis of higher order moments. Importantly, 
the joint PDF of local density, P(5,t;5',t'), evaluated at the same Lagrangian 
position but at the different times t and t' from the ellipsoidal collapse model 
exhibits a large amount of scatter. The mean relation between 5 and 5' does 
fail to match the one-to-one mapping obtained from spherical collapse model. 
Moreover, the joint PDF P(5; 6) from the ellipsoidal collapse model shows a 
similar stochastic feature, both of which are indeed consistent with the recent 
result from N-body simulations. Hence, the local approximation with ellipsoidal 
collapse model provides a simple but a more physical model than the spherical 
collapse model of cosmological PDFs, consistent with the leading-order results of 
exact perturbation theory. 
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Subject headings: cosmology: theory - galaxies: clustering - galaxies: dark matter 
- large-scale structure of universe - methods: analytical 

1. INTRODUCTION 

The large-scale structure of the universe is thought to be developed by the gravitational 
instability from the small Gaussian density fluctuations. In a universe dominated by dark 
matter, the evolution of mass distribution is entirely governed by gravitational dynamics. 
While luminous objects such as the galaxies and the clusters are subsequently formed by 
complicated processes including gas dynamics and radiative process, the dark matter distri- 
bution is the most fundamental product in the hierarchical clustering of the cosmic structure 
formation. In particular, the clustering feature of dark matter distribution is directly ob- 
served via weak gravitational lensing effect (Bartelmann & Schneider 2001 for review and 
references there in). Thus, the statistical study of the large-scale mass distribution provides 
a useful cosmological tool in probing the formation mechanism of dark matter halos, as well 
as the observed luminous distribution. 

In principle, all the statistical information of dark matter distribution is characterized 
by the probability distribution functions (PDFs) for mass density fluctuation and velocity 
fields, 5 and v. Among these, the one-point PDF for density field, P(5), is one of the 
most fundamental statistical quantities and because of its simplicity, there has been nu- 
merous theoretical studies on its evolution from a Gaussian initial condition. From the 
analytical study of one-point PDFs, Kofman et al. (1994) first calculated the PDF using the 
Zel'dovich approximation. For a perturbative construction of one-point PDF, Bernardeau & 
Kofman (1995) and Juszkiewicz et al. (1995) introduced the Edgeworth expansion to derive 
the analytic formula for PDFs. On the other hand, based on the tree-level approximation, 
Bernardeau (1992b) constructed the PDF from the generating function of the vertex func- 
tion. Interestingly, the vertex function in the tree-level approximation is obtained as an 
exact solution of the spherical collapse model. The effect of smoothing has been later incor- 
porated into his calculation and the analytic PDF was compared with N-body simulations 
(Bernardeau 1994b). Following these results, Fosalba & Gaztahaga (1998a) proposed to use 
the spherical collapse model for the prediction of higher-order moments beyond the tree-level 
approximation. The most remarkable point in their approximation is that the leading-order 
prediction exactly matches the one obtained from the rigorous perturbation theory and the 
correction for next-to-leading order is easily computed by solving the simple spherical col- 
lapse dynamics. Further, the spherical collapse approximation is recently extended to the 
prediction of one-point PDF (Scherrer & Gaztanaga 2001, see also Protogeros & Scherrer 
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1997). The approximation has been tested against N-body simulations and a good agreement 
was found even in a non-linear regime of the density perturbation. 

On the other hand, from the numerical study, Kayo, Taruya & Suto (2001) recently 
showed that the log-normal PDF quantitatively approximates the one-point PDF P(S) in the 
N-body simulations emerging from the Gaussian initial condition, irrespective of the shape 
of initial power spectra. The log-normal PDF has been long known as an empirical model 
characterizing the N-body simulations (e.g., Kofman et al. 1994; Bernardeau & Kofman 1995; 
Taylor & Watts 2000) or the observed galaxy distribution (e.g, Hamilton 1985; Bouchet et 
al. 1993; Kofman et al. 1994). Recently, a good agreement with the log-normal model has 
been reported in the numerical study of weak lensing statistics (Taruya et al. 2002a), and 
an attempt to clarify the origin of the log-normal behavior has also made (Taruya et al. 
2002b). Mathematically, the log- normal PDF is obtained from a one-to-one local mapping 
between the Gaussian and the non-linear density fields. Thus, the successful fit to the 
N-body simulation might be interpreted to imply that the evolution of local density can 
be well-approximated by the one-to-one local mapping. Indeed, the analytic PDF for the 
spherical collapse approximation can also be expressed as a one-to-one local mapping via the 
spherical collapse model. 

The above naive picture has been critically examined by Kayo, Taruya & Suto (2001) 
evaluating the joint PDF P(Si,ti, ^2,^2), i-e., joint probability of the local density fields Si 
and S2 at the same comoving position but at the different times ti and t%. It has been found 
that a large amount of scatter in the relation between Si and 82 shows up and their mean 
relation significantly deviates from the prediction from the one-to-one log-normal mapping. 
Although this might not be surprising, the good agreement between the log-normal and the 
simulated PDFs becomes more contrived and difficult to be explained in a straight forward 
manner. 

Definitely, the failure of the one-to-one local mapping somehow reflects the non-locality 
of the gravitational dynamics. That is, the evolution of local density cannot be determined 
by the initial value of the local density only. Rather, the local density can be expressed as 
multi-variate functions of local density and other local quantities such as velocity, gradient 
of local density, velocity-divergence and so on. Furthermore, initial values of these local 
quantities are randomly distributed over the three-dimensional space. As a consequence, 
even if the dynamics is deterministic, the relation between the evolved and the initial local 
density inevitably becomes stochastic. In a sense, the situation is quite similar to the non- 
linear stochastic biasing of galaxies, i.e., the statistical uncertainty between galaxies and 
dark matter arising from the hidden variables (e.g., Dekel & Lahav 1999; Taruya, Koyama 
& Soda 1999; Taruya & Suto 2000). Then, taking account of this stochastic nature, the 
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crucial issue is to construct a simple but physically plausible model of one-point PDF, at 
least, consistent with the N-body simulations in a qualitative manner. Further, the influence 
of stochasticity on the evolution of local quantities should be investigated. 

The purpose of this paper is to address these issues starting from a general theory of 
evolution of one-point PDF. In particular, we derive an evolution equation for the density 
and the velocity-divergence PDFs and consider how the stochastic nature of the local density 
field emerges. Due to the incompleteness of the equations, any theoretical approach using the 
evolution equations for PDFs generally becomes intractable. However, under the local ap- 
proximation that the evolution of density field (or velocity-divergence) is entirely determined 
by the local dynamics with finite degrees of freedom, we explicitly show that the analytical 
solution for the evolution equations is consistently constructed. Based on this approxima- 
tion, the one-point PDFs for the density and the velocity-divergence are computed from the 
ellipsoidal collapse model, as a natural extension of the one-to-one mapping of the spher- 
ical collapse approximation. Further, the stochastic nature arising from the multi-variate 
function of local quantities is explicitly shown evaluating the joint PDFs of the local density 
and/or the velocity-divergence. The influence of this effect is discussed in detail comparing 
with the spherical collapse approximation. 

This paper is organized as follows. In section 2, we consider a general framework to 
treat the evolution of one-point PDFs and derive the evolution equations for the Eulerian and 
the Lagrangian PDFs (Sec. 2. 2 and 2.3). Then, adopting the local approximation, consistent 
solutions for these equations are obtained (Sec. 2. 4). Further, the stochastic nature of the 
evolution of PDFs is quantified evaluating joint PDFs (Sec. 2. 5). As an application of these 
general considerations, in section 3, we give an explicit evaluation of one-point PDFs adopting 
the spherical and the ellipsoidal collapse models as representative Lagrangian local dynamics. 
The qualitative features of the results are compared with the perturbative analysis presented 
in appendix B or the previous N-body study. Finally, section 4 is devoted to the conclusion 
and the discussion. 

2. EVOLUTION EQUATION FOR PROBABILITY DISTRIBUTION 

FUNCTION 

2.1. Preliminaries 

Throughout the paper, we treat dark matter as a pressure-less and non-relativistic fluid. 
Assuming the homogeneous and isotropic background universe, the density field 5(x,t), the 
peculiar velocity field v(x,t) and the gravitational potential for the fluid follow the equation 
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of continuity, the Euler equation and the Poisson equation as follows: 

I + V«i + *W = o, (i) 

dv 1 1 

— + -(vV)v + Hv = — V0, (2) 
ot a a 

V 2 (f) = 47rGp m a 2 5, (3) 

where a is the scale factor of the universe, H(= a/a) denotes the expansion rate and p m is 
the background mass density. 

While we are mainly interested in characterizing the large-scale structure on the basis 
of statistical properties of the density field 5(x,t) and the velocity field v(x,t) as described 
in section 1, the dynamical evolution of such quantities is not solely determined locally. 
Hence, we must introduce the other local quantities characterizing the non-local properties 
of gravitational evolution, e.g., the gradient of local density V<5, the velocity tensor div/dxj, 
and so on. Let us denote these variables by 

dv ' \ 

8(x, t), v(x, t), V5(x, t), -^{x, t), ■ ■ ■) , (4) 

and define the joint PDF, P(^F;t), which gives a probability density that the quantity J- 
takes some values of (<5, v, V<5, • • •) at the time t. In terms of this, the one-point PDF of the 
density fluctuations P(5; t) is given as 

P(S;t)= ! 1]^^;*)- (5) 

and similarly the one-point PDF of the dimensionless velocity divergence, P(9;t), is 

P(0',t) = J Yl dFi P(F;t), (6) 



where 9 = V ■ v/(aH). 



In general, a statistical characterization of the large-scale structure is coordinate-dependent. 
Physically, there are at least two meaningful coordinates, i.e., the Lagrangian and the Eule- 
rian coordinates. While the Eulerian coordinate is fixed on a comoving frame, the Lagrangian 
coordinate is fixed on fluid particles and follows the motion of the fluid flow. Hence, as time 
goes on, high density regions in the Lagrangian space occupy larger volume than those in 
the Eulerian space. We thus consider both the Lagrangian PDF P L and the Eulerian PDF 
P E , defined in the Lagrangian and the Eulerian coordinates, q and x, respectively. The 
corresponding expectation values, (■ ■ -) L and (• • -) E are also introduced. 
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In the following subsection, we first consider the Lagrangian PDF and derive the evolu- 
tion equation. Then we derive the evolution equation for Eulerian PDF in section 2.3. The 
evolution equations derived here are not yet closed because of the unknown functions. In sec- 
tion 2.4 we discuss an approximate treatment using the local dynamics model, which enables 
us to obtain a closed form of evolution equation and to construct a consistent solution. 



2.2. Lagrangian one-point PDF 

To derive the evolution equation for the Lagrangian one-point PDF, we introduce an 
arbitrary function of local density, g(S), and consider the evolution of its expectation value, 
(g(S(q,t))) L evaluated in a Lagrangian frame. The differentiation of this expectation value 
with respect to time becomes 

|<^M))>, = | / UdF, g(6) P L (F;t) = J dSg(6) J^OM), (7) 

i 

since the explicit time dependence of the quantity (g(S)) L only appears through the La- 
grangian PDF. In the second equation, the integration is performed over the variable T 
except for the local density 5. On the other hand, the function g implicitly depends on time 
through the evolution of Lagrangian local density 5(q,t). Denoting the Lagrangian time 
derivative by d/dt = d/dt + v/a ■ V, the expectation value of the quantity dg/dt becomes 

The right-hand-side of the above equation can be expressed by integrating by part as 

i i 

P L (S;t)\. (9) 

8 ) 

Here, the quantity [A] B denotes the conditional mean of A for a given value of B: 

[A] B = I J] <LFiAP(F\B) (10) 

with the function P(J-\B) being the conditional PDF for a given B, i.e., P(J-\B) = 
P{F)/P{B). 





d5 


- 

OS \ 


dt 
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Now, recalling the fact that g(S) is an arbitrary function of local density 5, equation 
(7) must be equivalent to equation (8) in the Lagrange frame. The comparison between 
equation (7) and equation (9) then leads to the following evolution equation: 

'dS 



dt 



P L (S;t) 



0. 



11) 



Similarly, one can derive the evolution equation for the Lagrangian PDF of the velocity 
divergence P L (9;t): 

~dQ~ 



dt 



p L (e-,t) 



o. 



(12) 



2.3. Eulerian one-point PDF 

The evolution equation for the Eulerian one-point PDFs can also be derived by repeating 
the same procedure as above, but the resultant expressions are slightly different from those 
of the Lagrangian PDFs. The time derivative of the expectation value (g(S) ) E becomes 

d_ 

dt 



-(g(5(x,t))) E = J d5 g(5)^P E (5;t). (13) 



As for the expectation value of dg/dt, with a help of the Lagrangian time derivative, we 
obtain 

In the above expression, the first term in the second equation reduces to the same expression 
as in equation (9) just replacing the Lagrangian PDF with the Eulerian PDF. The second 
term in the second equation is further rewritten as 

= i (V •{».,(«(*,*)}>, - H(e 9 (5(x,t))), 

(X 



-H J Y[<LF i9 (6)eP E (r;t), (15) 



where we have used the fact that the first term in the second line vanishes because of the 
isotropy. Then, using the definition of the conditional mean (10), equation (14) can be 
summarized as follows: 

"d6 



P B (6;t) ) hH[e(t)],P,(6;t)}. (16) 



8 



Hence, the relation (dg/dt) E = d{g(S)) E /dt, not the equation (dg/dt) E = d(g(5)) E /dt, leads 
to the evolution equation for the Eulerian one-point PDF P(5,t). From equations (13) and 
(14), we obtain 



Comparing equation (17) with equation (11), the only difference between the Eulerian and 
the Lagrangian PDF is the term H[6]sP E in the right-hand-side of the above equation, which 
represents the change of the measure along the fluid trajectory in the Eulerian coordinate. 

Finally, the evolution equation of the Eulerian one-point PDF P E {9; t) is also derived in 
a similar way: 



Note that the zero-mean of the velocity divergence (9) E = is always guaranteed from 
equation (18), which is easily shown by integrating the above equation over 9 directly. 



The evolution equations in the previous subsection are rather general and in deriving 
these we did not specify the dynamics of fluid evolution. In this sense, they are not closed 
until functional forms of the conditional means [d8/dt] s , [9\ 5 and [dO/dt] g are specified. In 
other words, these equations require the additional equations governing the evolution of the 
conditional means. However, an attempt to obtain the closed set of evolution equations 
suffers from serious mathematical difficulty, so-called closure problem, which is well-known 
in the subject of fluid mechanics and/or plasma physics (e.g., Chen et al. 1989; Goto & 
Kraichnan 1993). Similar to the BBGKY hierarchy, if one derives the evolution equations 
for the conditional means, there appear new unknown conditional means. Thus, we must 
further repeat the derivation of evolution equation for the new quantities. Continuing this 
operation, one could finally obtain the infinite chain of the evolution equations, which is 
generally intractable. 

Instead of the exact analysis tackling the difficult problem, we rather focus on an approx- 
imate treatment of the evolution of one-point PDFs, where the solutions for the evolution 
equations are consistently constructed. To implement this, we adopt the local approximation 
that the evolution of the local density S and the velocity- divergence 9 is described by the 
Lagrangian dynamics with finite degrees of freedom, whose initial conditions are character- 
ized by the initial parameters, p = (pi,P2, • • -Pn), given at the same Lagrangian coordinate. 




(17) 




(18) 



2.4. Local approximation 



As will be explicitly shown in the next section, for instance, if the spherical collapse model 
is adopted as Lagrangian local dynamics, the evolution of local density is characterized by 
the single variable, which can be set as the linearly extrapolated density fluctuation, 5i. If 
adopting the ellipsoidal collapse model, the dynamical degrees of freedom reduce to three, 
representing the principal axes of initial homogeneous ellipsoid, Ai, A2 and A3. Thus, in 
this approximation, the density fluctuations can be expressed as 5 = f(p,t), and using this 
expression, the velocity-divergence is given by 6 = — (df/dt)/H(l + /) from the equation 
of continuity (1). Within the local approximation, provided the initial distribution function 
Pi(p), the form of the conditional means can be completely specified and it can be expressed 
in terms of the functions -Pj(p) and f(p,t). 

Let us first consider the Lagrangian PDF. In this case, the formal expressions for the 
conditional means [dS/dt] s and [d6/dt] e are given by 



~d5~ 




. dt . 


s P L (S;t)J 


'de' 




_dt_ 


e~ P L (0;t)j 



l[d,,, P I {p) d ^-5 f(p.t)). 



dt 



n^*<p)s{- 



df 



H(l + f) dt 



5d[0 + 



df 



H(l + /) dt 



(19) 



.(20) 



With these expressions, the evolution equations (11) and (12) become a closed form and the 
consistent solutions can be constructed as follows: 



P L (S;t) 

P L (9;t) 



I Udp.P^p) S D (S-f(p,t)), 
J i 



df(p,t) 



dt 



(21) 
(22) 



The proof that the above equations indeed satisfy the evolution equations (11) and (12) can 
be easily shown by differentiating equations (21) and (22) with respect to time. For the PDF 
of the local density, one has 

8 



8_ 

dt 



P L (S;t) 



l[dp t P I (p)j- t 5 D (5-f(p,t)) 

i 

Udp.P^p) 



dj_ 

dt 







85 



5 D (5-f(p,t)) 



from the property of the Dirac's delta function. In the above equation, the operator 8/ 85 
in the last line can be factored out and one can use the expression of conditional mean (19). 
Then, the time derivative of the one-point PDF P L (5;t) is rewritten as 
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which coincides with the evolution equation (11). As for the velocity-divergence PDF P L (9; t), 
we consistently recover the evolution equation (12) with a help of equation (20): 



\{dp l P I {p) 



d 



df\ d 



d_ 



dt 



dt [H(l + f)dt j de 
P L (0;t) 



>D 



e + 



df 



H(l + /) dt 



(24) 



The approximate solution of the Eulerian one-point PDFs are also obtained similarly, 
but the factor 1/(1 + 5) must be convolved with the Lagrangian PDF due to the presence of 
inertial term (r.h.s of eqs.[17][18]): 



P E (S;t) 
P E (0;t) 



1 + 5 



WdPiP^p) 5 D (5 - f(p,t)), 



1 df 
H{l + f(p,t)} dt 



P,t) 



(25) 
(26) 



Note, however, that these PDFs do not satisfy the following conditions: normalization con- 
dition (1) E = 1 and zero means (S) E = and (9) E = 0. This fact simply reflects that the 
conservation of Eulerian volume cannot be always guaranteed, in contrast to the conserva- 
tion of Lagrangian volume ensured by the mass conservation. As pointed out by Fosalba & 
Gaztanaga (1998a) (see also Protogeros & Scherrer 1997), we re-scale the relation between 
5 and f(p,t) as follows: 



6 = g(p,t) = N E {l + f(p,t)} -1 ; N E (t) = J \\d Pl - 



P;(P) 



+ f(p,ty 

Adopting this re-definition, the conditional means [dS/dt\ s and [d9/dt] e become 

ir*^/n^^)|(p.*)M*-.(P,o), 

( P ,t) 5 D (e-h( P ,t)), 



(27) 



d5_ 
~dl 

de 



dt 



where we define 



tt , PAp)_ dh 
P.(0;t)J U dPi 1+g dt 



h(p,t) 



dg 



H{\ +g)dt 



(28) 
(29) 

(30) 



Further, the conditional mean [9]s can be expressed as 



Ms 



1 \d8 



(31) 



H(l + 5) [dt 



from the equation of continuity (1). Then, the solutions of equations (17) and (18) becomes 



One can easily show that equations (32) and (33) satisfy the evolution equations (17) and 
(18), with the correct normalization and the zero mean. 

The above solutions for Eulerian PDF can be regarded as a generalization of the previous 
study based on the Zel'dovich approximation (Kofman et al. 1994) and/or the spherical 
collapse model (Scherrer & Gaztahaga 2001, see also Protogeros & Scherrer 1997). Note, 
however, that the general expression of velocity-divergence PDF P E (9;t) differs from the 
one obtained previously. While the factor 1/(1 + 5) is convolved in the integral in equation 
(33), the resultant expression by Kofman et al. (1994) obviously omitted this factor. In our 
prescription, the PDF P E (9;t) is constructed from the evolution equation, which basically 
relies on the equation of continuity. This means that, even for the velocity-divergence PDF, 
the factor 1/(1 + 5) is necessary to ensure the mass conservation. In fact, the expression (33) 
can also be obtained from the joint PDF P E (5, 9; t) integrating over the local density 5 (see 
eq.[36]). Although the velocity-divergence PDF in the previous study has been obtained in 
a rather phenomenological way, our present approach using the evolution equations might 
be helpful in constructing the consistent PDFs. 

Nevertheless, even at this point, the solutions of evolution equations should be regarded 
as formal expressions. In order to evaluate the PDFs explicitly, we need to specify the 
Lagrangian local dynamics. In other words, the quantitative prediction for PDFs using the 
local approximation crucially depends on the choice of the local dynamics. We will see in 
the section 3 that the explicit evaluation of P E (S; t) and P E {9; t) adopting the spherical and 
the ellipsoidal collapse models shows several noticeable differences. 



P E {S;t) 



f Y[d Pi P,(p) 5 D (5-g(p,t)) 



(32) 




(33) 



2.5. 



Joint PDF 



So far, we have restricted our attention to the one-point PDF for the single local variable, 
5 or 9. In our treatment of the local approximation, the expressions for PDFs are general 
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forms irrespective of the Lagrangian local dynamics. However, it should be emphasized that 
if the local dynamics is characterized by more than the two initial parameters, qualitative be- 
havior could be significantly changed from the local dynamics with single degree of freedom. 
The point is that the relation between initial parameters and the evolved result 5 or 9 cannot 
be described by a one-to-one mapping. Accordingly, the relation between 8 and 9 becomes 
no longer deterministic. Moreover, the failure of deterministic property also appears in the 
time evolution of such local variables. It is therefore important to discuss the stochastic 
nature of 5 and 9 arising from the dynamical evolution. To characterize this, we consider 
the joint PDF. Within the local approximation, one can construct a consistent solution of 
Eulerian joint PDF between 5 and 9 evaluated at the same time, P E (5,9;t). Further, the 
Lagrangian joint PDF for the density field evaluated at the same Lagrangian position but 
at the different times, P L (5,t; S',t') can also be obtained. 

The evolution equation of P E (5,9;t) can be derived through the expectation value of 
an arbitrary function g(S,9). Repeating the same procedure as described in section 2.3, we 
obtain 



l(g(6(x,t),e(x,t)))> 
^g(5(x,t),9(x,t)) 



d8d9 g(S,e)^P E (6,9;t), 



/ 



= / dSd9 g(5,9) 



d_ 



d5 
~dl 



P E (S,9;t) 



8,6 




d9 
~dl 



P E (5,9;t) + H9P E (5,9;t) 



5,9 



Then, these two equations lead to the evolution equation for P E (S,9;t): 



^P E (6, 9; t) - H9^- s {(1 + 5)P E (5, 0;t)} + ^ 



d9_ 
~dl 



P E (S,9;t) =H9P E (5,9;t), (34) 



5,8 



where we used the relation [dS/dt]s,e = —H0(1 + 5). 

The construction of the consistent solution in the local approximation is almost parallel 
to the case of the Eulerian one-point PDF in section 2.4. The formal expression of the 
conditional mean [d9/dt]s t e is 



d9 
~dl 



5,6 



P^WlTTS /II** f M*- g (P«->«M)), (35) 
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and the corresponding solution of equation (34) becomes 

1 



P B (5,6;t) = 



1 + 5 



H d Pt P 7 (p) 5 D (5 - g(p, t))5 D {9 - h(p, t)). 



(36) 



Note again that the above solution exactly recovers the one-point PDFs P E (S; t) and P E (9; t) 
integrating equation (36) over 9 and 5, respectively (see eqs.[32] [33]). 

The evolution equation for Lagrangian joint PDF is also obtained by the time derivative 
of the expectation value of the arbitrary function, g(5(q,t),5(q,t')) as described in section 
2.2: 



d_ 
~dt 



(g(5(q,t),5(q,t>))) 



g(5{q,t),S(q,1?)) 



d5d5'g(5,5')^ t P L (5,t;5>,t>), 

d ( , ( . xf ,,,d5(q,t) 
g(5(q,t),5(q,t)) 



d5(q,t) 

I 



d6d5'g(5,6')^- s 





'd6~ 




dt_ 



dt 

P L (S,t; 5',t') 



S,6> 



The evolution equation of Lagrangian joint PDF is 

d 



dt 



P L (5,t;5',t') + 





~d5~ 


- 

05 { 


_dt_ 



P L (5,t;S',t') =0. 



(37) 



5,8' 



The conditional mean in the local approximation is expressed as 



d5 
di 



6,8' 



Pl(6. 



J 11^ P +P) 6 D (6-f(p,t))5 D (6>-f(p,t>)). 



Recalling that the joint PDF satisfying the evolution equation (37) should be invariant under 
the transformation, (S,t) <-> (8', If), the solution consistent with the boundary condition 
P L (5, t'; 5', t') = P L (5; t')8 D (8 - 5') becomes 



P,(8,t;8',t') = JUdptP^SoiS-fip^SniS'-f&t')). 



(38) 



Notice that if the local Lagrangian dynamics is described by a single parameter, the integral 
over the initial parameter pi in equation (38) can be formally performed. The resultant 
expression includes Dirac's delta function, leading to the one-to-one mapping between 5 and 
5'. On the other hand, in cases with the multivariate initial parameters, one cannot generally 
perform the above integral and the Dirac's delta function is not factored out, leading to the 
stochastic nature of local density fields. 
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One might further consider the evolution of Eulerian joint PDF P E (8,t;S',t'), which 
has been indeed shown in N-body simulations by Kayo et al. (2001). The derivation of 
evolution equation itself is an easy task, but, the formal solution in the local approximation 
suffers from difficulties due to the presence of advective term, which might be related to an 
important effect on the non-local nature of fluid flows. Since even the Lagrangian joint PDF 
shows several important features, one can expect that the qualitatively similar behavior to 
the N-body results can be seen from the Lagrangian joint PDF. Hence, we will postpone to 
analyze the Eulerian joint PDF P E (5, t; 5', t') and instead focus on the Lagrangian joint PDF 
P L (5,t';5',t>). 



3. DEMONSTRATION AND RESULTS 

Now we are in a position to give an explicit evaluation of the one-point PDF based on 
the local approximation discussed in the previous section. For an illustrative purpose, we 
adopt the spherical and the ellipsoidal collapse models as simple and intuitive Lagrangian 
local dynamics. After briefly describing the basic equations of these models in section 3.1, 
we compute the Eulerian one-point PDFs P E (S) and P E {9) and discuss the qualitative differ- 
ences arising from the choice of Lagrangian local dynamics in section 3.2. In particular the 
stochasticity in the multi-variate function of local density or velocity-divergence is examined 
in detail by evaluating the joint PDFs, P L (5, t; 5', t') and P E (S; 6) from the ellipsoidal collapse 
model. 



3.1. Models for Lagrangian local dynamics 

First consider the simplest case of the Lagrangian local dynamics, in which the evolution 
of local quantities is determined by the mass inside a sphere of radius R collapsing via self- 
gravity: 

d 2 R GM „ 4vr 



(Jf , R2 M = — pR i = const, (39) 

where M is the mass inside the radius and p represents the local density. This spherical 
collapse model can be re-expressed as the evolution equation for density fluctuations 5, 
given by 5 = ~p/p m — 1 = (a/R) 3 — 1 as follows: 

d 2 5 nTT d5 4 1 fd5\ 2 3 TT2f ^ „ ... 

with the quantity f2 m being the density parameter, f2 m = 8nG p m /(3H 2 ). As Fosalba & 
Gaztanaga (1998a) stated, this equation is regarded as a shear-less and irrotational approx- 
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imation of fluid equations, since one can derive the following equation from equations (1) to 
(3): 

§ + < - iih (f ) 2 = ^ + * + - ^) • 

with a help of the Lagrangian time derivative. Here the quantities and lo^ respectively 
denote the shear and the rotation given by 

°« = + ( 42 ) 



1 1 


f dvi 


+ ^ 


2aH 1 






1 i 


( dvi 


_ 3^ 


2a# 1 







Wy " 2^fe~^J- (43) 

In the spherical collapse approximation (40), the density fluctuations S depend only on a 
single initial parameter Si, i.e., the linear fluctuation at an initial time, if the decaying mode 
of the linear perturbation is neglected. Note that the spatial distribution of the initial density 
field is randomly given and thereby the parameter Si is regarded as a random variable. We 
assume that the initial parameter Si obeys a Gaussian distribution: 

P t (Si) = -jL- e -(*/ ff ') a / 2 , (44) 

V27T Gl 

where the variable 07 means the rms fluctuation of the linear density field Si, i.e., of = (Sf). 

The ellipsoidal collapse approximation which is another Lagrangian local dynamics, 
describes the evolution of the uniform density ellipsoid. In contrast to the spherical collapse 
model, the evolution of S is governed by the dynamics of the half length of principal axes 
ctj (i = 1,2,3) characterizing the density ellipsoid. According to Bond & Myers (1996), we 
have 

d 2 „ fl + S h 



l t i ai = -^Gp m at y— h —6 + \ cxt ,iJ , (45) 

f°° dr 2 

bi = a x a 2 a z I 7 - r - — — - -, (46) 

Jo («i i+^)IL 1/2 3 



and the relation between S and a,- becomes 



a 3 



S = 1. (47) 

Here, the variable A extj j quantifies the external tidal effect, which is required for the consis- 
tency with the Zel'dovich approximation in a linear regime (Bond & Myers 1996): 

. D(t) ( Aj — ^ 1 ~*~ ^ 2 ^~ ^ 3 ) ; linear external tide , , , 

Acxt.i — \ V 3 J (48) 

| bi ; nonlinear external tide , 
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where D(t) is the linear growth rate. The variables A« represent the initial parameters of 
principal axes, and in terms of these, the initial conditions reduce to 

<*i(t ) = a{to){l-D(t )Xi}, (49) 

^Ui(to) = a(t ){l - D(t )Xi} - a(t )D(t )Xi. (50) 
at 

In contrast to the spherical collapse model, one can regard this model as an approximation 
of the fluid equations taking account of the influence of shear but neglecting the rotation: 

d 2 5 „„d5 4 1 fd5\ 2 2 „ /3 



■W + 2H Tt-- 3 —An) =^(1 + *)^* + ^! («> 

<Tij = — 3 dij. 

SH \ a>i oti a 2 a 3 J 

Note that similar to the spherical collapse approximation, the three initial parameters \ 
are regarded as the random variables. When the initial condition of density field is assumed 
to be a Gaussian random distribution, the expression for the initial parameter distribution 
P/(Aj) is analytically obtained as follows (e.g., Doroshkevich 1970; Bardeen et al. 1986): 

3375 / T 2 T \ 

P^i) = -^rs exp -3-4, + 15-^ (Ax - A 2 )(A 2 - A 3 )(A 1 - A 3 ), (52) 

where the quantities 1\ and l 2 denote 1\ = \\ + A 2 + A3 and 1 2 = A1A2 + A 2 A 3 + A 3 Ai, 
respectively. 

Based on these Lagrangian local models, we numerically calculate the PDFs assuming 
the Einstein-de Sitter universe (f2 m = 1, = 0), in which the linear growth rate D is 
simply proportional to the scale factor a. For a better understanding of the later analysis, 
in Figure 1, we plot the evolution of local density 5 from the ellipsoidal collapse model for 
some initial conditions (e,p) given by e = (Ai — \ 3 )/(25i) and p = (Ai + A 3 — 2A 2 )/(2<5j). The 
results are then depicted as a function of linearly extrapolated density 5i = Ai + A 2 + A 3 and 
are compared with the one from the one-to-one mapping of spherical collapse model (solid). 
Figure 1 shows that the local density of the ellipsoidal collapse model generally takes a larger 
value than that of the spherical collapse model. Further, the variety of evolved density for 
fixed 81 suggests that a large amount of scatter will appear in the joint PDF P L (5,t;5',t') 
and the resultant one-point PDFs P E (S) and P E (6) cannot, in general, coincide with those 
obtained from the spherical collapse model. 
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3.2. Results 

In computing the PDFs from the above local collapse models, one may practically en- 
counter the case when the local density infinitely diverges at finite elapsed time for some 
regions in the initial parameter space, which has not been treated in previous section. To 
avoid the unphysical divergences, we must restrict the integral in the PDFs to the initial pa- 
rameter space V(t), in which the local density 5 does not diverge. Indeed, this modification 
slightly affects the normalization condition for both the Lagrangian and the Eulerian PDFs, 
i.e., ( 1 ) LiE 7^ 1. Although this does not alter any qualitative features of PDFs, we consider 
some modifications to keep the correct normalization and adopting this procedure in ap- 
pendix A, and the results for one-point and joint PDFs are presented below. Note, however, 
that the perturbation calculation discussed in 3.2.1 is free from the serious divergences and 
within the perturbative treatment, one can rigorously develop the local approximation for 
PDFs. 



3.2.1. One-point PDFs 

Let us show the results of the one-point PDF. Figure 2 plots the one-point Eulerian 
PDFs of the local density (top) and the velocity-divergence (bottom) evaluated at various 
linear fluctuation values o\. In both panels, the thick lines represent the results obtained from 
the ellipsoidal collapse model with linear external tide, while the thin lines denote the PDFs 
from the spherical collapse models. In computing the PDFs, the Lagrangian local dynamics 
are numerically solved with the various initial conditions p in the initial parameter space 
V. Then, weighting by the PDF of the initial parameter Pi(p), the PDFs are constructed 
by binning the evolved results of the density S and the velocity-divergence 0, together with 
appropriate convolution factors (see Appendix A). 

As expected, the overall behaviors of both PDFs in Figure 2 are qualitatively similar, 
irrespective of the Lagrangian local models. As increasing the linear fluctuation value 07, 
while the density PDF P E (S) extends over the high-density region 5 ^> 1, the velocity- 
divergence PDF P E (0) is negatively skewed and it extends over the negative region <C — 1. 
In looking at the differences in each local model, we readily observe several remarkable 
features. First, the density PDFs computed from both the spherical and the ellipsoidal 
collapse models almost agree with each other. At first glance, this seems to contradict with 
a naive expectation from the local dynamics in Figure 1. However, one might rather suspect 
that the agreement in density PDFs is accidental, due to the distribution of initial parameters 
Pf(Aj), which is, at least, consistent with a naive picture that joint PDF P L (5, t; 5', t') from the 
ellipsoidal collapse model exhibits a large mount scatter and the mean relation between 5 and 
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5' significantly deviates from one-to-one mapping of spherical collapse model (see Sec. 3. 2. 2 
and Fig. 4). On the other hand, the velocity-divergence PDFs from the ellipsoidal collapse 
model exhibit longer non-Gaussian tails, compared with those obtained from the spherical 
collapse model. The deviation between both models in PDF P E (9) becomes significant as 
increasing the value 07. Interestingly, in the non-linear regime o\ > 1, tails of PDF P E (0) 
from the spherical collapse model show the opposite tendency, i.e., the amplitude decreases 
as increasing 07. 

In order to characterize the qualitative behaviors more explicitly, we perturbatively 
solve the evolution equations for both the spherical and the ellipsoidal collapse models. The 
differences are then quantified evaluating the higher order moments of one-point statistics 
for the local density and velocity-divergence. In appendix B, based on the formal solution of 
PDFs in section 2.4, perturbative calculations of local collapse models are briefly summarized 
and the solutions up to the fifth order are presented. The resultant expressions for the higher 
order moments of density and velocity-divergence are obtained as a series expansion of linear 
variance of, up to the two- loop order for the variance and the one-loop order for the skewness 
and the kurtosis: 

a 2 = (5 2 ) s = of + s 2 , 4 of + s 2 , 6 <rf + •••, (53) 

53 = ^ = S 3 ,o + S 3 , 2 of + •••, (54) 

54 = - — ^-r. = £4,0 + 5*4,2 of H (55) 



for the local density and 



a 2 = (6 2 ) E = a 2 + 4 A al + 4 fi af + ---, (56) 
T 3 = ^l£ = T3,o + T 3 , 20 f + •••, (57) 



T 4 = (g4) % 34 =T 4 , + T 4 , 2 of + ■■■ (58) 

for the velocity-divergence. Then, all the coefficients in the above expansions yield the 
rigorous fractional number and Table 1 displays a summary of the results. The calculation 
in spherical collapse model essentially reproduces the non-smoothing results obtained by 
Fosalba & Gaztanaga (1998a,b). Note, however, that the discrepancy has appeared in the 
higher order correction of velocity-divergence (c.f., eq.[12] with 7 = of Fosalba & Gaztanaga 
1998b). Perhaps, in computing the velocity-divergence moments, Fosalba & Gaztanaga 
(1998b) incorrectly used the cumulant expansion formula for S listed in Fosalba & Gaztanaga 
(1998a). Further, we suspect that they erroneously replaced the convolution factor 1/(1 + 5) 
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in Eulerian expectation value with 1/(1 + 0). On the other hand, in our calculation, moments 
(9 n ) E are rigorously computed according to the definition (B6), with a help of the velocity- 
divergence PDF (33) with equation (30). Hence, the present calculation is at least consistent 
with the local approximation in section 2.4 and we believe that no serious error has appeared 
in present result. 

Apart from this discrepancy, one finds that the leading-order (tree-level) calculation of 
skewness S^o and kurtosis S^q in both models exactly coincides with each other. While 
the differences in the higher order correction for local density are basically small, consistent 
with Figure 2, the results in velocity-divergence exhibit a large difference, especially in the 
variance Oq. Figure 3 summarizes the departure from the leading-order perturbations for the 
variance (top), the skewness (middle) and the kurtosis( bottom), each of which is normalized 
by the leading term. Clearly, the higher order corrections for variance show the significant 
difference between the spherical and the ellipsoidal collapse models, although the model de- 
pendence of the external tide in ellipsoidal collapse is rather small. Remarkably, the one-loop 
correction is negative in the spherical collapse model and thereby the quantity ag does 
not monotonically increase. This behavior indeed matches with the non-monotonic behavior 
of velocity-divergence PDF seen in Figure 2. In this sense, the perturbation results success- 
fully explain the numerical results of PDF. This fact further indicates that in a Gaussian 
initial condition, the influence of non-sphericity or effect of shear could be negligible in the 
one-point statistics of local density, while it alters the shape of the PDF P E (9), which might 
be a natural outcome of the multivariate local approximation. 



3.2.2. Joint PDFs 

Next we focus on the joint PDFs. Left panel of Figure 4 shows the Lagrangian joint 
PDF P L (5(z = 0);5(z = 9)) from the ellipsoidal collapse model, evaluated at the present 
time z = with various linear variances 07. On the other hand, right panel of Figure 4 
represents the results fixing the linear fluctuation value to u\ = 2 at present, but at different 
output times. Although Figure 4 does not rigorously correspond to the N-body results 
obtained by Kayo et al. (2001) (c.f. Fig.7 in their paper), the qualitatively similar features 
can be drawn in several points. First, the scatter between 5(z = 0) and 5(z = 9) becomes 
broader as increasing the time and the linear variance (top to bottom in left panel). Second, 
the nonlinearity between the initial and the evolved density indicated from the curvature 
of the conditional mean [5 (z)] s( z =9) (solid) also tends to increase as time elapses (top to 
bottom in right panel). The one-to-one mapping obtained from the spherical collapse model 
(short- dashed) is very different from the mean relation [S(z)]s( z =9), but their mean relations 
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basically reflect the qualitative behavior of local dynamics in Figure 1. That is, the evolved 
results of local density in the ellipsoidal collapse tends to take a larger value than that in 
the spherical collapse. Moreover, recall the fact that both the initial and the final PDFs of 
local density P E (S) show a good agreement between the spherical and the ellipsoidal collapse 
model (see Fig. 2). This is indeed the same situation as in the N-body simulation; apart from 
the detailed differences, a simple model of PDFs provides an essential ingredient for the 
stochastic nature of N-body results. In this sense, the local approximation with ellipsoidal 
collapse models can be regarded as a consistent and physical model of one-point statistics, 
which successfully explains the N-body simulations. 

Finally, using the ellipsoidal collapse model with linear external tide, we examine the 
Eulerian joint PDF of local density and velocity-divergence evaluated at the same time, i.e., 
P E (8; 9). In Figure 5, contour plots of joint PDF P E (8; 9) for various linear fluctuation values 
oi are depicted as function of —9 and 5. This is the so-called density-velocity relation, which 
might be of observational interest in measuring the density parameter Q m from the velocity- 
density comparison through the POTENT analysis (e.g., Bertschinger & Dekel 1989). Along 
the line of this, theoretical works based on the Eulerian perturbation theory have attracted 
much attention, as well as the N-body study (e.g., Bernardeau 1992a; Chodorowski & Lokas 
1997; Bernardeau et al. 1999). Based on the solution of the local approximation (36), one can 
easily calculate the perturbation series of velocity-density relation [9]s as function of local 
density and density- velocity relation [S]e as function of velocity-divergence, the leading-order 
results of which are expected to coincide with previous early works in the non-smoothing 
case, exactly. Beyond the perturbation analysis, Figure 5 reveals the general trend of the 
stochastic nature in the velocity-density relation. As increasing 07, the scatter becomes much 
broader and the conditional means [8]g (dot- dashed) and [9}s{solid) does not coincide with 
each other. Of course, the one-to-one mapping obtained from the spherical collapse model 
(short- dashed) fails to match the both conditional means. These qualitative behavior is in 
fact consistent with the N-body results by Bernardeau et al. (1999) and the present model 
provides a simple way to derive the non-linear and stochastic velocity-density relation. 



4. CONCLUSION AND DISCUSSION 

In this paper, starting from a general theory of evolution of one-point PDFs, we derived 
the evolution equations for PDF and within a local approximation, consistent formal solu- 
tions of PDF are constructed in both the Lagrangian and the Eulerian frames (see eqs. [21] [22] 
for Lagrangian PDFs and eqs. [32] [33] with eqs. [27] [30] for Eulerian PDFs). In order to reveal 
the stochastic nature arising from the multivariate Lagrangian dynamics, we further con- 



- 21 - 



sider the Eulerian joint PDF P E (8, 9; t) (eq.[36]) and the Lagrangian joint PDF P L (5, t; 5', t') 
(eq.[38]). Then, adopting the spherical and the ellipsoidal collapse models as representative 
Lagrangian local dynamics, we explicitly evaluate the Eulerian PDFs, P E (S) and P E (0), as 
well as the joint PDFs. The results from the ellipsoidal collapse model show several dis- 
tinct properties. While the PDF P E {5) almost coincides with the one-to-one mapping of 
the spherical collapse model, the tails of velocity-divergence PDF P E {6) largely deviate from 
those obtained from the spherical model. These behaviors have also been confirmed from 
the perturbative analysis of higher order moments. On the other hand, evaluating the La- 
grangian joint PDF of local density, P L (5, t; 5', t'), a large scatter in the relation between the 
initial and the evolved density fields was found and their mean relation fails to match the 
one-to-one mapping of spherical collapse model. This remarkably reproduces the same situ- 
ation in the N-body simulation. Therefore, the local approximation with ellipsoidal collapse 
model provides a simple and physically reasonable model of one-point statistics, consistent 
with the leading-order results of exact perturbation theory. 

Since the present formalism described in section 2 is quite general, the approach does 
not restrict its applicability to the pressure-less cosmological fluid. Rather, one may apply 
to the various fluid systems in presence or absence of gravity. As mentioned in section 
2.4, however, the applicability or the validity of local approximation of PDFs, in principle, 
sensitively depends on the choice of Lagrangian local dynamics. In the last section, simple 
and intuitive examples were examined for the illustrative purposes. The results indicate that 
the multivariate Lagrangian dynamics rather than the local model with a single variable can 
describe various statistical features of fluid evolution including the stochastic nature. 

Perhaps, a straightforward extension of the present treatment is to include the effect 
of redshift-space distortion or projection effect, which is practically important for proper 
comparison with observation. Before addressing this issue, however, remember the most 
primarily importance of the smoothing effect. While the appropriate prescription for top-hat 
smoothing filter does exist in the local approximation with the spherical collapse model (e.g., 
Bernardeau 1994b; Protogeros & Scherrer 1997; Fosalba & Gaztanaga 1998a), the smoothing 
effect on the approximation using ellipsoidal collapse model still needs to be investigated. 
This step is in particular a crucial task in order to construct a more physical prescription 
for one-point statistics of cosmic fields and the analysis is now in progress. The results will 
be presented elsewhere (Ohta, Kayo & Taruya, in preparation). 

We are grateful to Y.Suto for reading the manuscript and comments. I.K acknowledges 
the support by Takenaka-Ikueikai Fellowship. This work is supported in part by the grand- 
in-aid for Scientific Research of Japan Society for Promotion of Science (No. 1470157). 
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A. ON NORMALIZATION CONDITION OF PDFs 

In computing the PDFs from the local approximation with the spherical or the ellip- 
soidal collapse model, one may practically encounter the case when the evolved density at a 
finite elapsed time infinitely diverges for some regions in initial parameter space. To avoid 
the unphysical divergences, we restrict the integral over the entire initial parameter space 
appearing in the PDFs to some regions V(t), in which the local density 5 does not diverge. 
This modification slightly alters the normalization condition for both the Lagrangian and 
the Eulerian PDFs in section 2.4 and 2.5, which should be corrected in the following re- 
normalization procedure. 

First of all, the initial distribution function P/(p) should be re-normalized as 

P,(P) = , Pi(p) . (Al) 



Mt) V 



Then, the modification of Lagrangian PDFs is to replace the initial distribution P/(p) with 
P/(p), together with the change of the integral region. For example, the one-pint PDF P L (S) 
in equation (21) is modified as follows: 

P L (S;t)= [ l[dp t P I (p)5 D (5-f(p,t)). (A2) 

Mt) i 

On the other hand, for the Eulerian PDFs, notice the fact that the re- normalization (Al) 
also affects the relation between 5 and f(p,t) (see eq. [27]), which must be modified as 

l+g(p,t) = N E (t){l + f(p,t)} ; N E (t)= f II d * i I'tf + \ ■ ( A3 ) 
Using this relation, the renormalized one-point PDFs (32) and (33) respectively become 

^;*) = TX7 / UdPi Pi(p) S D (6 - g{p,t)) (A4) 
for the PDF P E (5) and 

p ^L$ d ^A e+ ^l +k ) (A5) 

for the PDF P E {6). Here the variable k is given by 

k -wJ rit M d -mrwf t p - iP) ' (A6) 
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under the approximation that the time evolution of V is neglected, V/V <C f/(l + /). 
Similarly, the renormalized Eulerian joint PDF P(S, 9; t) corresponding to the expression 
(36) becomes 

* «> = TTs l m n *■ p - W M* - » fo (" + 77(1^7) I + *) ' (A7) 



B. CUMULANTS FROM SPHERICAL AND ELLIPSOIDAL COLLAPSE 

MODELS 

In this appendix, based on the local approximation with spherical and ellipsoidal col- 
lapse models, we briefly summarize the essence of the perturbation analysis for higher order 
moments of local density and velocity-divergence. The details of the calculation procedure 
including the effect of smoothing will be presented elsewhere (Ohta, Kayo & Taruya, in 
preparation). Here, we only present the results in non-smoothing case. 

Let us write down the evolution equation for ellipsoidal collapse model. In an Einstein- 
de Sitter universe, equation (45) becomes 

2 d 2 ai adai 3 / 1 + 8 .h \ . . 

° drf ~ 2 da = -2* (— + 2 6 + ) ■ (B1) 

In this case, the linear growth rate D is proportional to the scale factor a. Thus, the half 
length of principal axis ctj can be expanded by a power series of a: 




Note that the initial conditions (49) and (50) state 7^ = \. Hence, substituting the 
expression (B2) into (Bl) and solving the evolution equation perturbatively, the coefficient 
7^ is systematically determined order by order and is expressed as the function of \. Thus, 
the perturbative expansion of local density 5 given by (47) is also expressed in terms of Aj : 

f(\ t ,a) = - 1 = V 6®(\i) a\ (B3) 

with the corresponding boundary condition being 5^ = Ai + A 2 + A 3 . 

Provided the function /(Aj, a), one can calculate the n-th order moments of 5 and 6 as 
well as the normalization factor N E) according to the one-point PDFs of local approximation, 
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(32) and (33): 



N E = 



JV £ {l+/(Ai)} 

r% = j de r Pm / !, . , -J-^ _ ^^)"P,(A ( )(B6) 

where the function P I (Xi) denotes the PDF of initial parameter Aj given by (52). 

Below, we separately present the perturbation results in each model. For the results of 
higher order moments, i.e., variance, skewness and kurtosis given by (53)-(58), numerical 
values of the perturbation coefficient are summarized in table 1 and the departure from the 
leading-order results are depicted in figure 3. 



B.l. Ellipsoidal collapse model with linear external tide 

If adopting the ellipsoidal collapse model with linear external tide (see eq. [48]), the 
perturbative expansion of local density, 5^ up to the fifth order becomes 

5 (1) = S h (B7) 

6(2) = 2T^ + 2T Jl ' (B8) 

567 ' 4725 ' 405 ' V ; 

. U) 55805 r4 952144 „ 345088 . T 12368 t2 , r>in , 

S {4) = 5? H 5 2 Ji H Si J 2 H Jh (BIO) 

130977 ' 2480625 ' 16372125 363825 1 y 1 

c(5) 213662 r5 237342074 r3 93363344 r2 

5 (5) = 5? H 5, 3 Ji H <J T Jo 

729729 ' 621928125 ' 3192564375 ' 

52865818 2 135052 

+ 638512875^ + 34827975 JlJ2 ' (BU) 

where 5/ denotes the linear fluctuation density given by 5i = Ai + A 2 + A 3 . Here we introduced 
the quantities J\ = x 2 + xy + y 2 and J 2 = (x — y) (2x + y) (x + 2y) with the variables x and 

y being x = Ai — A 2 and y = A 2 — A3, respectively. Substituting the above expansion results 
into (B4), the normalization factor of Eulerian PDF becomes 

69668 , 

N ° = 1 + 3898125*' ? + "" < B12 > 
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Further, using the equations (B5) and (B6), one can obtain the variance, the skewness and 
the kurtosis defined in (53)-(58). For the variances of local density and velocity-divergence, 
the perturbative correction up to the two- loop order 0(af) becomes 



57137 

33075 ( 
8747 



469828713881 

111739753125 
154583563 



-ai —————ay + ..., (B13) 

= af + —a?+ < + ■■■■ (B14) 

1 11025 ' 165540375 1 K J 

As for the skewness and the kurtosis, we obtain the results up to the one-loop order 0(af ): 

34 646404856 



S 3 



+ 



63669375 



< + 



60712 210688932175742 
+ 



1323 



782178271875 



-a + 



for the local density and 



T 3 
T 4 



26 333940984 



7 

12088 
441 



63669375 



< + 



+ 



17145801103334 
156435654375 



-a + 



(B15) 
(B16) 

(B17) 
(B18) 



for the velocity-divergence. 



B.2. Ellipsoidal collapse model with non-linear external tide 

In the case of the ellipsoidal collapse model with nonlinear external tide, the perturbative 
solution 5^ up to the fifth order becomes 

= 5 h (B19) 

6(2) = 2T^ + 2T Jl ' (B20) 

. m 341 ^ 338, T 92 /rtn ,. 

S (3 > = Sf H 5,Ji H J 2 , (B21) 

567 945 3969 

_ m 55805 r4 485288 r2 234088 . 429728 t2 /T10oA 

<5 (4) = 5, H 5? Ji H 5,J 2 H J?, (B22) 

130977 ' 1091475 ' 4584195 10696455 1 K J 

213662 - 292398464 rS 64182728 r2 T 

S {5 > = Sf H <J?Ji H 5? J 2 

729729 ' 638512875 ' 893918025 ' 

6541246 e t2 828974992 .„ , 

+ 5;Ji + Ji J 2 - B23 

59594535 1 96364363095 K ' 

Then the normalization factor is 
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The variances of local density and velocity-divergence up to the two-loop order become 



439 

245 
145 



3143785639 

695269575 
1708470649 



= cr, 



-a 



147 1 1158782625 
The results for the skewness and the kurtosis can be reduced to 

34 1041064 



for the local density and 



S3 

T 3 
T 4 



101871 



■<7f + 
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for the velocity-divergence. 
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B.3. Spherical collapse model 

In the case of spherical collapse model, the perturbative solutions of local density f(5i) 

just corresponds to those obtained from the ellipsoidal collapse model simply setting J± — 
and J2 = 0, since the three initial parameters of principal axis are identical. The cumulants 

are then calculated in similar way to the previous subsection except for the initial parameter 
PDF (44). The resultant normalization factor is 

N E = 1 + —a] + -^-at + ■■■. (B31) 

21 ' 43659 ' V ; 

The variances of 5 and 9 become 

1909 4 344439415 6 

a — Ui H <7i H <7 1 +•••, rJoz) 

1 1323 1 107270163 1 ' y ! 

a 2 e = a}-— <jf+™*!®-af + .... (B33) 

e 1 U7 1 1324323 1 y ' 

The skewness and the kurtosis respectively becomes 

34 1026488 9 .„ , 

5 3 = — + a? + • • • , B34 

3 7 101871 ' ' V ; 

60712 22336534498 2 

54 " "1323" + 83432349 + " ' ' (B35) 
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for the local density and 



T 3 
T 4 



26 613936 



7 101871 1 
12088 10081115810 9 



441 



83432349 



(B36) 
(B37) 



for the velocity-divergence. 
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Table 1. Coefficients of perturbative correction 
for the density and the velocity-divergence in 
one-point statistics from the spherical and the 
ellipsoidal collapse models (see eqs.[53]-[58]). 





SCM t 


ECM1 t 1 


ECM2 * 2 


S2,4 
«2,6 


1.44 
3.21 


1.73 
4.20 


1.79 
4.52 


-S'3,0 
5*3,2 


4.86 
10.08 


4.86 
10.15 


4.86 
10.22 


S'4,0 
5*4,2 


45.89 
267.72 


45.89 
269.36 


45.89 
270.79 


„e 

b 2,4 

JS 

6 2,6 


-0.075 
0.24 


0.79 
0.93 


0.99 
1.47 


^3,0 
^3,2 


-3.71 
-6.03 


-3.71 
-5.24 


-3.71 
-5.30 


^4,0 
^4,2 


27.41 
120.83 


27.41 
109.60 


27.41 
110.62 



t Spherical collapse model 
^ Ellipsoidal collapse model with linear external tide 
Ellipsoidal collapse model with non-linear external tide 



31 




Fig. 1. — Evolution of local density from the ellipsoidal collapse model for various initial 
parameters, (e, p), where e is ellipticity and p is prolateness, given by e = (Ai — A3) / (25i) and 
p = (Ai + A 3 — 2A 2 ) / (25 1), respectively. The evolved results are then plotted as a function of 
linearly extrapolated density, 5i = Ai + A2 + A3. For comparison, we also plot the one-to-one 
local mapping from the spherical collapse model (solid). 
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Fig. 2. — Eulerian one-point PDF, P E (S)(top) and P E (9) (bottom) for the ellipsoidal collapse 
model with linear external tide (thick lines) and for the spherical collapse model (thin lines). 
In each panel, solid, long-dashed and short-dashed lines represent the results in cases with 
the linear rms fluctuation 07 = 0.5, 1 and 2, respectively. 
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Fig. 3. — Departures of the variance, the skewness and the kurtosis of local density (left) and 
the velocity-divergence(ng/i£) from the tree-level perturbation results. While the variances a 2 
and erg are calculated up to the sixth order of the linear rms fluctuation 0/, the perturbation 
results for the skewness S3, T 3 and the kurtosis £4, T 4 are obtained up to the one-loop order, 
0(af) (see Appendix B in detail). In each panel, the dashed and the dotted lines indicate the 
results obtained from the ellipsoidal collapse model with linear external tide and nonlinear 
external tide, respectively. For comparison, the results from the spherical collapse model are 
also shown in solid lines. 
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(5(z = 9) (5(z = 9) 

Fig. 4. — Contour plot of Lagrangian joint PDF P L (8(z),8(z = 9)) obtained from the 
ellipsoidal collapse model with linear external tide. The left panel shows the results fixing 
the redshift to z — 0, while the right panel represents the evolution of joint PDF fixing the 
linear variance to o\ = 2. In each panel, solid lines indicate the conditional mean [5(2)]$( z= g) 
computed from joint PDF. For comparison, the one-to-one local mapping obtained from 
the spherical collapse model are also plotted in short-dashed lines. Left: o\ = 0.5 (top), 
a i = l(middle) and o\ = 2(bottom). Right: z = 2(top), z = l(middle) and z = O(bottom). 
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Fig. 5. — Contour plot of Eulerian joint PDF P E (S;6) as functions of —9 and 5 from the 
ellipsoidal collapse model with linear external tide. In each panel, the solid and the dot- 
dashed lines respectively denote the conditional means —[9]s and [5]g. The short-dashed 
lines represent the relation obtained from the spherical collapse model: o\ = 0.5(top); o\ = 
l(middle); oi = 2(bottom). 



